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Abstract 

In this paper, we give an analytical treatment to study the behavior of the col- 
lapse and the revival of the Rabi oscillations in the Jaynes-Cummings model (JCM). 
The JCM is an exactly soluble quantum mechanical model, which describes the in- 
teraction between a two-level atom and a single cavity mode of the electromagnetic 
field. If we prepare the atom in the ground state and the cavity mode in a coherent 
state initially, the JCM causes the collapse and the revival of the Rabi oscillations 
many times in a complicated pattern in its time-evolution. In this phenomenon, the 
atomic population inversion is described with an intractable infinite series. (When 
the electromagnetic field is resonant with the atom, the nth term of this infinite se- 
ries is given by a trigonometric function for y/nt, where t is a variable of the time.) 
According to Klimov and Chumakov's method, using the Abel-Plana formula, we 
rewrite this infinite series as a sum of two integrals. We examine the physical mean- 
ings of these two integrals and find that the first one represents the initial collapse 
(the semi-classical limit) and the second one represents the revival (the quantum 
correction) in the JCM. Furthermore, we evaluate the first and second-order per- 
turbations for the time-evolution of the JCM with an initial thermal coherent state 
for the cavity mode at low temperature, and write down their correction terms as 
sums of integrals by making use of the Abel-Plana formula. 

1 Introduction 

The Jaynes-Cummings model (JCM), which described the interaction between a two- 
level atom and a single electromagnetic field mode, was originally proposed for examining 
spontaneous emission in 1960s [Ij. This model is derived from applying the rotating wave 
approximation to an electric dipole coupling. In the interaction term of the Hamiltonian of 
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the JCM, the photon creation operator accompanies the atomic de-excitation operator and 
the photon annihilation operator accompanies the atomic excitation operator. Because 
the JCM is an exactly soluble quantum mechanical model, it is investigated theoretically 
by researchers in the field of quantum optics eagerly [21 El EJ- 

If we initially prepare the atom in the ground state and the cavity mode in a coherent 
state, the JCM causes the collapse and the revival of the Rabi oscillations many times 
in a complicated pattern in its time-evolution and this phenomenon is regarded as the 
evidence of the quantum nature of the electromagnetic field E] • (This phenomenon was 
confirmed experimentally in 1980s [7].) Thus, the demonstration of the collapse and the 
revival of the Rabi oscillations in the JCM gives the foundation to Planck's thought [8]. 
That is, the collapse and the revival of the Rabi oscillations in the JCM tells us that the 
photon's energy is equal to hi/, where h represents the Planck's constant and v represents 
the frequency of the photon, so that excitation of the photons shows discreteness. 

Recently, the JCM has been studied from a new viewpoint by the researchers in the 
field of quantum information science. The JCM is often used for investigating the evolu- 
tion of entanglement between the atom and the single mode cavity field [91 [10] . The lower 
bound of entanglement between the two-level atom and the thermal photons in the JCM 
is also discussed [TTj . The JCM can be applied to the realization of quantum computa- 
tion [12J. So-called sudden death effect (disappearance of entanglement of two isolated 
Jaynes-Cummings atoms in a finite time) is predicted [T31 [EU H5] , and it is experimentally 
demonstrated [16]. Thus, some researchers in the field of quantum information science 
think that the JCM has to be studied from the new viewpoint. 

When we discuss the JCM, we often have to handle an intractable infinite series. 
(If the electromagnetic field is resonant with the atom, the nth term of this infinite 
series is given by a trigonometric function for y/nt, where t represents a variable of the 
time.) For example, the atomic population inversion in the collapse and the revival of 
the Rabi oscillations in the JCM is described by this infinite series. In the thermal JCM, 
whose initial state of the cavity mode is given by a thermal equilibrium state, the atomic 
population inversion is written as a similar intractable infinite series, as well. 

In Ref. [17] . Klimov and Chumakov discuss the thermal JCM and evaluate the atomic 
population inversion, which is described by the intractable infinite series. They change 
this intractable infinite series into a sum of two integrals by making use of the Abel-Plana 
formula mentioned in Ref. [18] . By carrying out numerical calculations, they find that the 
first integral represents a semi-classical limit (the initial collapse) and the second integral 
represents a quantum correction (quasi-chaotic behavior). 

In this paper, we give an analytical treatment to study the behavior of the collapse 
and the revival of the Rabi oscillations in the JCM, according to the method proposed by 
Klimov and Chumakov. For applying the Abel-Plana formula to the infinite series that 
describes the atomic population inversion, we replace an inverse of a factorial 1/n! with 
an inverse of the gamma function l/r(n + 1) and perform the analytical continuation 
on the complex plane as 1/T(z + 1). (This prescription is a new key point of this paper 
as compared with Ref. [1TJ .) After giving this step, using the Abel-Plana formula, we 
write down the atomic population inversion as a sum of two integrals. We examine the 
physical meaning of these two integrals and find that the first integral represents the 
initial collapse (the semi-classical limit) and the second integral represents the revival 



2 



(the quantum correction) in the JCM. In this paper, we clarify that we can separate the 
quantum correction from the semi-classical limit in the solution of the JCM. 

Furthermore, we evaluate the first and second-order perturbations for the 
time-evolution of the JCM, whose initial state of the cavity mode is given by a thermal 
coherent state at low temperature [T9l [20] . [The expansion parameter for the perturbation 
theory is given by 9(13) ~ exp(— (3e/2), where j3 = (k^T)' 1 , e = his, j3e 3> 1 and v is the 
frequency of the cavity field. A rigorous definition of 9(0) is given by Refs. [19] and [20].] 
We obtain the first and second-order correction terms of the atomic population inversion 
and rewrite them as sums of integrals, using Klimov and Chumakov's method and the 
Abel-Plana formula. 

After deriving integral forms of the atomic population inversion for both resonant 
and off-resonant cases and their thermal perturbation corrections, we examine effects of 
detuning and low temperature against the collapse and the revival of the Rabi oscillations 
numerically. 

This paper is organized as follows. In section [2j we give a review of the JCM and 
its collapse and revival of the Rabi oscillations. In section [3j we apply Klimov and 
Chumakov's method to the atomic population inversion, which is represented by the 
intractable infinite series. (In this section, we consider the zero-temperature case.) We 
rewrite it as a sum of two integrals with the Abel-Plana formula. In section HJ we consider 
the time-evolution of the JCM at low temperature. Preparing the cavity field in a thermal 
coherent state initially, we evaluate the first and second-order perturbation corrections of 
the atomic population inversion. Using the Abel-Plana formula, we rewrite them as sums 
of integrals. In section [51 we examine properties of the integrals obtained in sections [3] 
and H] in detail by numerical calculations and analytical methods. In section [HI we give a 
brief discussion. In Appendix [AJ we give a derivation of the Abel-Plana formula, which 
plays an important role in Klimov and Chumakov's method. In Appendix [Bj we give 
some remarks about precise techniques for carrying out the numerical calculations shown 
in section [5j 



2 The JCM and its collapse and revival of the Rabi 
oscillations 

The JCM is a quantum system, which consists of the single two-level atom and the single 
mode of the electromagnetic field. Its Hamiltonian is given by 

h 

H = —ujq(T z + hucva + hn(a + a + cr_a^), (1) 

where % = h/2n, a± = (l/2)(a x ±ia y ) and [a, at] = 1. The Pauli matrices (a iy i = x,y,z) 
are operators of the atom, and a and a< are an annihilation and a creation operators of 
photons, respectively. Moreover, we assume the coupling constant k to be real. 
We can divide H given by Eq. (OQ) as follows: 

H = h(Ci + C 2 ), 

d = u[(l/2)a g + a!a], 

C 2 = n(a + a + a_a f ) - (Au/2)a z , (2) 
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where Au = u — luq. Because [Ci,C2] = and we can diagonalize C\ at ease, we take 
the following interaction picture. We write a state vector of the whole system in the 
Schrodinger picture as \ips(t)). We define a state vector in the interaction picture as 
\ipi(t)) = exp(iCit)|$3(i)). [We assume |^i(0)) = |^s(0)).] The time-evolution of |-0i(^)> 
is given by |^(t)) = U(t)\ik(0)), where U(t) = exp(-iC 2 t). 

We give the basis vectors for the state of the atom and the photons of the cavity field 
as follows. (In this section, we consider only the zero-temperature case.) We write the 
ground and excited states of the atom as two-component vectors, 



\9)a = " > I*>a = n , ( 3 ) 





where we assume that |<?)a and |c)a are eigenvectors of a z with eigenvalues —1 and 1, 
respectively. (The index A stands for the atom.) We describe the number states of the 
photons as |n) P (n = 0, 1, 2, ...), which are eigenstates of a^a. (The index P stands for the 
photons.) 

Describing the atom's Pauli operators by 2 x 2 matrices, we can write down U (t) as 
follows: 



where 



/ o x i * sin(v / D + kH) 



sin(V-Di) 

uqi = ik -j= — a, 

\/D 



sin(VD + kH) , 

■ -n I 



Urn = IK- 



r— % . sin(\/I)i) 
«n = cos{VDt) + -Aw — ^_ y , (5) 

D = (Aw/2) 2 + reVa. (6) 

Here, we define the initial state of the whole system as follows. We assume that the 
atom is in the ground state \g)\ at t = 0. Moreover, we assume that the cavity mode is 
in the coherent state \a)p at t — 0, where 

00 a n 

\a) P = exp(-|a| /2) ^ -=\n) P , (7) 
n=o v n! 

and a is an arbitrary complex number. (We provide that 0! = 1.) From now on, to let 
the discussion be simple, we assume a to be real. Writing the initial state of the whole 
system as \ipi(0)) = | c?)a|«)p, the time-evolution \tpi(t)) is given as follows: 

t*oo «oi 1 / \ _ ( u 01 \a) P 



and 
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Figure 1: A graph of (a z (t)) (t G [0, 20ir}) for k — 1, Au = and a = 4. In the numerical 
calculations of (a z (t)) defined in Eq. (fill) , the summation of the index n is carried out 
up to n = 100. Looking at this graph, we notice that the initial collapse time is order of 
unity and the period of the revival is approximately equal to 8tt. 



Thus, the probability that we observe \g)A at the time t is given by 

P 9 (t) = \A(gMt))\ 2 



exp(— a 2 ) — r[ cos ~ 



oo a 2n 



n=0 



nl 



; v /(Aw/2) 2 + nK 2 t) 



+ (Aw/2) 2 



sin 2 (J(Aw/2) 2 + nK 2 t)_ 



(Aw/2) 2 + m 2 1 
The atomic population inversion is given by 

(<j z {t)) = Tr P (V>i(0k^i(*)) = 1 - 2P g {t). 



(9) 



(10) 



Especially, in the case where Au = 0, that is, the electromagnetic field is resonant with 
the atom so that the energy gap is equal to the frequency of the photons, we can replace 
\K,\t with t by letting the time t be in units of and we obtain 



oo a 2n 



(v z (t)) = -exp(-a 2 ) V — cos(2y/nt). 



11) 



We plot (a z (t)) (t e [0, 20vr]) for k = 1, Aa; = and a = 4 in Fig. UJ 

Looking at Fig. [H we can observe the collapse and the revival of the Rabi oscillations 

clearly. In general, the larger |a| is, the more distinctly we observe the revival of the Rabi 

oscillations. From Fig. [H we can suppose that the time scale of the initial collapse is order 

of unity and the period of the revival is around 27r|a|. 

The time scale of the initial collapse and the period of the revival are explained as 

follows [2TJ [22] . Let us evaluate (u z (t)) defined in Eq. (FlT!) with assuming a 2 3> 1. 
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Writing the index of the summation as n = a 2 + 5n, because of the property of the 
Poisson distribution, the major contribution to the right-hand side of Eq. All]) comes 
from the terms with \5n\ < a 2 , so that we can neglect the terms with \5n\ > a 2 . We 
rewrite y/n as 



- a 2 + n 

n ~ — : — — . 12) 

2\a\ K J 

Then, we can obtain an approximate form (a z (t)) in Eq. (fflil as 

(a z (t)) ~ _exp(-a 2 )f;^^[exp(2z^^t)+exp(-2z^^t)] 

= — exp[a 2 (cos - — : — 1)1 cos(|a|£ + a 2 sin - — r ). (13) 

|a| ' \a\ 

In the right-hand side of Eq. (Tl3|) . exp[a 2 (cos(t/|a|) — 1)] represents an amplitude enve- 
lope of the wave, which causes a kind of the beat, and cos[|a|i+a 2 sin(t/|a|)] represents the 
Rabi oscillations. If we assume 0<£/|a|<Cl,we obtain cos(t/|a|) ~ 1 — (t 2 /2a 2 ), and the 
factor of the amplitude envelope approximates to exp[a 2 (cos(t/|a|) — 1)] ~ exp(— 1 2 /2), so 
that we understand that the initial collapse time is approximately equal to unity. More- 
over, the factor of the amplitude envelope exp[a 2 (cos(£/|a:|) — 1)] shows that the period 
of the revival is given by 27r|a|. 

Thus, Eq. (fT3"|) represents the collapse time and the period of the revival of the genuine 
(cr z (t)) given by Eq. f fTTj) well. However, to investigate properties of (o~ z (t)) precisely, 
changing the infinite series that appears in Eqs. (Q and (flOj) into a simple form is favorable. 
This is the motivation of this paper. 



3 Representing the atomic population inversion as a 
sum of two integrals for the zero-temperature case 

To rewrite the atomic population inversion (cr x (t)) given by Eqs. (jUJ) and ffTUj) . we use the 
following formula. If n\ and n 2 are integers and <f>(z) is a function which is analytical and 
bounded for all complex values of z such that ri\ < Re(z) < n 2 , then 

^(m) + 4>(ni + 1) + <f>{m + 2) + ... + (j)(n 2 - 1) + ^0(n 2 ) 

rn2 roc \ 

x [<p(n 2 + iy) - (f>(ni + iy) - <p(n 2 - iy) + 0(m - iy)]dy. (14) 

(This formula is described in Ref. [18J as an example. The author of Ref. [18] mentions 
only a suggestion to prove this formula. We give details of the derivation of this formula 
in Appendix E) Moreover, if we assume <fi(z) — > as Re(z) — > +oo, we obtain 

\m + ± m = r + i r ^ iy) ~ y t { : iy) d y . ^ 

ii 1 6 -L 
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This equation is called the Abel-Plana formula. 

The Abel-Plana formula given by Eq. f|T5|) changes an infinite series into a sum of two 
integrals. However, we cannot apply this formula to (cr z (t)) given by Eqs. and ffTUl) in 
a straightforward manner. Thus, we try to extend this formula in the following way. 

First, we consider a series, 

OO 1 

£rr/(" + <0, (16) 



n=0 n - 



where c is a real constant and the function f(x) is infinitely differentiable and bounded 
for any real value of x. Thus, we can rewrite f(n + c) as the Taylor series and we obtain 



oo -r oo -r oo m 

L ~> - » 1 » c c(m), 



E zr/(" + c) = E rr E ^r/ (m) H- (17) 



„=0 n! "0 n! m! 



[In Eq. (TTZD, /( m )(n) denotes (d m /dx m )f{x 

Second, using the property of the gamma function T(n + 1) = n\, we rewrite Eq. ((1 

as 

OO OO £,771 OO 



E ^/(» + o) = E b E rv^T/^ (")> 

where 



oo 

z-1 -t 



r(z) = / t^e-'dt. (19) 
Jo 

The gamma function T(z) given by Eq. (JT9l) converges absolutely only for Re(z) > 0. 
However, by analytical continuation, we can let T(z) be analytical everywhere in the 
complex plane except at z — 0, —1, —2, .... Because there are no points at which T(z) is 
equal to zero, 1/T(z) is analytical at all finite points of the complex plane. 

Third, we assume that f^ m \z)/Y{z + 1) is analytical and bounded for all complex 
values of z such that < Re(z) < +oo and f^ m \z)/T(z + 1) -> as Re(z) -» +oo for 
m = 0,1,2,.... Then, we can apply the Abel-Plana formula given by Eq. ( TT5l) to the 
right-hand side of Eq. ( 1T8|) and we obtain 

oo i oo m 1 roo f("i)( r ) 

EV^ + C ) = E^T^ (m) (0)+/ l-Mrdx 
„r,n\ ±^ n m! L 2 Jo r [X + 1 



77=0 '"■ 777 =0 

roo 



1 / (m) fa) f( m \-iy) y 
+% I- e**v-l [ r(l + iy) F(l -iy) l V 



2 JK ' Jo r(x + i) 



-2 f-i— ^{ Z^tM }^. (20) 

Jo e 2 ^-l x r(l + i?/) J y v ; 

Here, applying the formula given by Eq. f l20|) to Eq. (jUJ), we can change P g {t) repre- 
sented as the infinite series into a sum of integrals. At first, we define 

c=(Au/2k) 2 . (21) 

Then, defining f(z + c) as 

f(z + c) = \a\ 2{z+c) [cos 2 (^/IT7:\ n\t) + -^-sm 2 (^/z~T^\K\t)\, (22) 
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we can confirm that f(z + c)/T(z + 1) is analytical and bounded for all complex plane. 
Moreover, f(z + c)/T(z + l) ^ as Re(z) +00. [In the limit of Re(z) +00, |r(z + l)| 
increases more rapidly than any exponential functions of z.] Thus, we obtain 



00 \ \2{n+c) 



_^ |a;| 2(nH ■ , , 

P g (t) = exp(-a 2 )|a| 2c ^2 L - L ~ ] {cos 2 (y/n + c\n\t) H sin 2 ( v / nTc|/c| 

n=o n - n + c 

= exp(-a 2 )[(l/2) + / 1 (0) (t)-24 0) (t)], (23) 

where 

poo \cx\ 2x C 

I{ } (t) = / 1 1 [cos 2 {Vx + c + l\n\t) + -sin 2 (Vx + c + l\n\t)]dx. 



T(x +1) x + c + Z 



4°(f) = r-J-rM 



e-»-r--r(rTM [cos( ^ + i+i!/|K|i) 



c 



+ —7—- sin^^c + / + iy\K\t)]}dy. (24) 
c + I + ly v 

In the above equation, P g (t) is represented as a sum of two integrals, J °° dx and J °° dy 
[namely, /{ 0) (t) and 4°\t)}. Thus, using Eq. ( ITOl) . we can describe (cr z {t)) as a sum of 
these integrals, 

(a z (t)} = 1 - exp(-« 2 )[l + 2l{°\t) - 4/ 2 (0) (t)]. (25) 
In the case where Aw = 0, k = 1 and c = 0, we obtain 

P s (t) = exp(-a 2 )[- + — - - cos 2 (y/xt)dx 

( \2iy 



roc 1 \a\ y 1 — 

\ +exp(- a 2 )[i + I|~ F ^_ T ^co S (2V5t)<fa 



-ri^ im{ f^) c ° s ( 2 ^)>*]. ( 26 > 

where we use 



-5^ -Im{ ^ ' , J rfy = -- + -exp(a 2 ). (27) 



2 Jo T(x+1) Jo e 2 *y-l l T(l + iy) J * 4 2 
[We can derive Eq. ( !27l) from Eq. ( ITBI) at ease.] Thus, we obtain 

(a a (t)> = -(1/2) exp(-a 2 ) + J x {t) + J 2 (t), (28) 

where 

yoo I 1 2a; 

J i{t) = — exp(— a ) / — cos(2y/xt)dx, 

Jo T(x + 1) 

yoo 1 \a\ 2iy I — 

J 2 (t) = 2eM-* 2 )J o i ^3 T Im{ f J r ^- y cos(2 V / ^t)}rf 2/ . (29) 
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4 The time-evolution of the JCM with an initial 
thermal coherent state 



4.1 The definition of the thermal coherent state 

In Refs. [19] and [20], a thermal coherent state is defined as an extension of the zero- 
temperature coherent state according to the thermo field dynamics. In this section, 
preparing the atom in the ground state and the cavity mode in the thermal coherent 
state initially, and letting the whole system evolve in time with the JCM, we calculate the 
atomic population inversion up to the second-order perturbation correction. The obtained 
first and second-order correction terms are represented as the intractable infinite series. 
In subsection 14. 4[ we rewrite these correction terms as sums of integrals by making use of 
the Abel-Plana formula. 

In the thermo field dynamics, we have to handle a space that is a direct product of 
the ordinary zero-temperature Hilbert space H and a so-called tilde space %. Thus, every 
number state of the photons \n) G H has a corresponding state |n) G "H, so that the 
orthogonal basis for the whole space in thermo field dynamics is described as {\n) <g> \m) G 
H <8) H : n, m G {0, 1, 2, ...}}. In the next paragraph, we give a definition of the thermal 
coherent state. 

First of all, we define creation and annihilation operators acting on 7-L and Iri as follows: 

[a, a'] = [a, a'] = 1, [a, a] = [a, a'] = 0. (30) 

We pay attention to the fact that the operators on % (namely a and a)) and the operators 
on % (namely a and a)) commute with each other. Next, we introduce the temperature 
by the following unitary transformation: 

Up = exp(-iG), (31) 

G = -z6 l (/3)(a t a t - aa), (32) 

cosh6(/3) = [1 - exp(-/3e)]" 1/2 , 

smh6{/3) = [exp(/3e) - I]" 1 / 2 , (33) 

where (5 = (k^T)' 1 and e = hu. We note that 0{(5) is real. To emphasize that G and 
Up are operators acting on both "H and H, we put an accent (a hat) on them. Then, 
annihilation operators are transformed as follows: 

a — > a{(5) = UpaUt, 
a —7- a{(5) = UpaUl. 

Writing down the zero-temperature vacuum as |0) <g) |0) 
vacuum is given by 

\0(/3)) = Up\0,0). 

Then, a thermal coherent state is defined as follows: 

\a, 7; 13) = exp[<W(/3) + 7*e^(/3) — a*a(/3) 



(34) 

= |0, 0) G H ® U, the thermal 

(35) 

-ia(P)]W)h ( 36 ) 
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where a and 7 are arbitrary complex numbers. From now on, for the simplicity, we assume 
a and 7 to be real. 

Writing the initial state of the whole system as IV'i(O)) = |<7)a|<*> T, P)p, the time- 
evolution of \ipi(t)) is described as 

\ uii\a,j,p)p J 

where u m and u\\ are defined in Eqs. (J5J) and (jSJ). Thus, the probability that we observe 
\g)h at the time t is given by 

^(0;*) = \u n \a,f,P)\ 2 - (38) 

To calculate P g ((3;t) by the perturbation theory, we rewrite \a, 7;/?) given by Eqs. fl34l) . 
f )35|) and f )36|) as follows: 

|a,7;/3) = exp[t/a(aa^ + 70^ — aa — ja)Up\Up\0, 0) 

= f/g exp[<W + 70^ — aa — 7a] |0, 0) 

= Up exp[a(a^ — a)] exp[7(a^ — a}\0, 0) 

= ^|«)|7>, (39) 

where \a) and (7) represent the (zero-temperature) coherent states in H and Iri, respec- 
tively. Thus, we can change Eq. fl38l) into the following form: 

P g ((3;t) = (amfyul&nUplaW). (40) 

Moreover, using Eqs. (jSJ) and ([6]), we can rewrite the Hermitian operator u^un as 

+ 9/1 , !~- x sin (\K\tvaJa + c) , , 

wli«n = cos 2 ( ac \t\/(tfa + c) + c iL 4 (41) 

a* a + c 



where c is defined in Eq. (12 ip . 

To compute P g (f3;t) approximately, we formulate the perturbation theory as follows. 
At first, we consider an arbitrary function g(x), which can be represented as a Taylor 
series about x — 0, 

00 

9(x) = E 9 (n) x n , (42) 

n=0 

where —00 < x < +00 and g( n ' is an arbitrary complex number for all n. [In Eq. (|42|) . 

g{n) is 

equal to (d n /dx n )g(x)\ x=0 . This notation is slightly different from f( m \n) given by 
Eq. (jTTJ) .] Next, we evaluate U^g(a)a+c)Up by the second-order perturbation theory at low 
temperature. We assume that the temperature T is quite low, so that j3e = e(k-oT)^ 1 3> 1. 
Then, from (3e ^> 1 and Eq. (I33l) . we understand < 9(f3) <C 1. Decomposing tanh#(/3) 
into the power series in the small parameter 9(0), 

tanh0(/9) = e~^ /2 

= e(f3)-(l/3)9(f3) 3 + 0(9 5 ), (43) 
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we obtain 9((3) ~ e _/3e//2 with neglecting terms of order 0(9 3 ). From now on, we consider 
the second-order perturbation theory with the parameter 9(0). 

From the above discussions, we can compute Upg(a)a + c)Up approximately as follows: 

Utg(ala + c)Up 

oo 

= g(a)a + c) + 9((3) £ # (n) [a t a t - aa, (a f a + c) n ] 



n=0 

oo 



+ (l/2)9((5) 2 ^ (n) [a t a t - aa, [atat - aa, (a f a + c) 1 



n=0 



+0(9 S ). (44) 

(In the above derivation, we use the Baker-Hausdorff theorem [2].) In the following 
subsections, we calculate the first and second-order perturbation corrections, respectively. 

4.2 The first-order correction 

In this subsection, we calculate the first-order correction given by Eq. (j4"4"j) . To obtain the 
commutation relation of (a)a) — aa) and (a^a + c) n , we use the following notation: 

A = a) a) — aa, B = a) a + c, C — a) a) + aa. (45) 

Then, we obtain the following relations: 

[A,B] = -C, [C,B] = -A. (46) 

Using Eq. (JIB"]) , we obtain 

[A,B] = -C, 
[A, B 2 ] = A-2BC, 
[A,B 3 ] = -C + 3BA - 3B 2 C, 

[i,S 4 ] = A - ABC + 6B 2 A - AB 3 C, (47) 

and so on. 

Here, we introduce the following operators: 

fx = era', v = aa. (48) 

Because 

A — pu — v, 6 = p, + v, (49) 
we can rewrite Eq. (j4"Tj) as follows: 



[A,B] 


= (B- 


i)A — « 


'l + B)u- 


BA, 


[A, B 2 ] 


= (B- 


i) 2 A - 


(l + B) 2 P 


- B 2 A, 


[A, B 3 ] 


= (B- 


i) 3 A- 


(l + Bfu 


- B 3 A, 


[A, B*] 


= (B~ 


iffi - 


(l + B) 4 u 


— B A A. 



(50) 
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Thus, we obtain 

[A,B n ] = (B-l) n fi-{B+l) n D-B n A for n = 1,2,3,.... (51) 
From Eqs. (150 p and (15 ip . we obtain the first-order correction as 

oo 

e{(3)Y,9 [n) [ci ] ^ -aa,(a^a + c) n } 

n=0 

= 6(PMB-l)fi-g(B + l)u-g(B)A). (52) 
4.3 The second-order correction 

In this subsection, we calculate the second-order correction given by Eq. ( I44"j) . From 
Eq. f[5"Tj) . we obtain 

[A, [A, B n ]] = R n + S n for n = 1,2,3,.., (53) 

where 

R n = [A,(B-l) n }(i-[A,(B + l) n )v-[A,B n )A, 

S n = (B-l) n [A,fx] + (B + l) n [A,u]. (54) 

Thus, we can divide the second-order correction up into the following two parts: 

-6{Pf £ g^R n + -6(P) 2 £ 9 {n) S n . (55) 

1 n=0 ^ n=0 

First, we evaluate the part including {S n } in Eq. fl55|) . Because of Eqs. ( |4"8"j) and fj4~9]) . 
we obtain 

[i,£] = [i, £] = -£>, (56) 

where 

D = a) a + a)a + 1. (57) 

[We note that the above D is different from .D defined in Eq. (jSJ).] Thus, we can write 
down S n as 

S n = -((B-l) n + (B + l) n )D. (58) 
Hence, we can rewrite the term that includes {S n } in Eq. (I5"5"j) as 

-±0(^ 2 (<7(B-1) + <7(B + 1))D. (59) 

Second, we evaluate the part including {i? ra } in Eq. ( 155]) . From Eq. (J5i|) . we can 
confirm i?o — -Ri — at ease. After slightly tough calculations, we obtain 

R 2 = -2[A,B]C, 

R 3 = -3[A,B 2 ]C + 3[A,B]A, 

R 4 = -4[A,B 3 ]C + 6[A,B 2 }A-A[A,B]C, (60) 
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and so on. Thus, we can write down the part that includes {R n } in Eq. ( 155]) as 

-1 oo 



n=0 



+ J2(n + 2)(n + 3)g^[A,B n+1 ]A 

1 n=0 
-r oo 

-J2(n + 2)(n + 3)(n + A)g {n+4) [A, B n+1 ]C + 



6 



n=0 



^(/3) 2 (E 



n=0 



A/! 



i \n oo ^ 

-J-F n fi - -Av 



n=0 



nl 



-\e{P) 2 (g{B - l)jM - g(B + 1)0 - g{B)A)A, 



where 



dx' 1 



x=B-l 



dx" 1 



x=B+l 



d n 

v - d^ 9{x] 



A. 



x=B 



(61) 



(62) 



[In the above derivation, we use Eq. ( 15T|) in an effective manner. The form of F n in 
Eq. ([62]) reflects Eq. (151}.] 

Here, we introduce the following operators: 



00 ( — "H™ A n 



n=0 



n! <ix r ' 



(63) 



Using the symbols e ±d//da: , we can rewrite the first term of Eq. (IBTj) as 

^(/3) 2 ([e" d/d ^(5 - - e- d/dx g(B + l)z> - e~ d/dx 'g{B)A]fi 
-\e d/dx g(B - l)/t - e d,dx g{B + 1)0 - e d,dx g(B)A}0 
Then, we apply the following technique to Eq. ((6 



e ±d l dx g{X) = £ 



n=0 



;-i) n d n 

n\ dx* 



x=X 



g{x±i) 



(64) 



(65) 



where X is an arbitrary operator. Thus, we can rewrite Eq. (1M]) as 

^(/3) 2 ([<?(£ - 2)/* - g(B)0 - g(B - l)A]fi 
-[g(B)fi - g(B + 2)0 - g(B + 1)A}0). 
Finally, from Eqs. fl55|) . fl59|) . f l6Tj) and fl66|) . we obtain the second-order correction as 

oo 

(l/2)6>(/5) 2 g (n) [a j ^ - ad, [a t S t - ad, (a f a + c) n ]] 



(66) 



n=0 



(l/2)fl(/3) 2 [s(£? - 2)^ 2 + g(B - 1)(-D - fiA - Ap) 

+g(B)(A 2 -Oft- fiO) + g(B + !)(—£> + 0A + AO) + g(B + 2)z> 2 ]. (67) 
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4.4 The atomic population inversion up to the second-order 
correction 

In this subsection, we compute P g ({3;t) denned in Eq. (|4"0l up to the second-order per- 
turbation. According to Eq. ( )4T1) . we define g(x) as follows: 

g(x + c) = cos 2 (\n\t\/x + c) H sin 2 (|/dtA/x + c). (68) 

X + c 

From Eqs. Q3DJ, (gU), ((52]) and (ETJ), we obtain P g (P;t) as 

P fl (/9;t) = P g (t) +9((3)PM(M + (l/2)#(/3) 2 Pf (/3;t) + 0(£ 3 ), (69) 

where 

PW(M = (a\m 9 (B-l)fi-g(B + l)u-g(B)A]\am, 
Pf\M = {a\mg{B-2)f? + g{B-l){-b-fiA-Af,) 

+g(B)(A 2 -vfi- (jlv) + g(B + 1)(-D + vA + AO) 

+g{B + 2)v 2 }\am. (70) 

First, we compute P^\(3;t). Remembering that \a) and I7) are coherent states on H 
and H respectively, we can write down pW(/3;t) as follows: 

P^\M = H<7l££exp( ' ' ' 



n=om=o 2 Vn\\"m\ 

x ^[^(a^a + c — 1) — g(a)a + c)]aW 
— [g(a^a + c + 1) — ^(a^a + c)]aa) \n)\fh) 
exp[-(a 2 + 7 2 )] 

00 oo 2n+l ~2m+l 

x(EE „, 7 ■ [g(n + c)-g(n + C +l)] 



00 00 ^2n— 1 ^2m— 1 

££<»-l)'<m-l)l b(n + c) - s(B + C - 1)1 



00 a 2n 



2a7exp(-a 2 ) £ — r[g(n + c) - #(n + c+ 1)] 

n=0 

207^*) -gW(*)], (7i) 



where 



00 a 2n c 

= exp(-a 2 ) V — -[cos 2 (a/c + n + /kit) H sin 2 (\/c + n + 1\k\ 

= exp(-a 2 )|«|- 2 ( c+ ^ Y \ a \ 2n+C+l [cos 2 (Vc + n + l\K\t) 
■sm 2 (Vc + n + l\K\t)} 



c + n + I 

exp(-a 2 ) -[cos 2 (|/«|tv / c+7) H -sm 2 (|/«|tv / e+7)] 

— C — (~ I 

+/ 1 (,) (t)-2/? ) (t)l for/ = 0,1,2,..., (72) 
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and l[ l \t) and 4°(t) are defined in Eq. flM}. We note that Q^(t) = P g (t). 

Second, we compute Pjj 2 \(3;i). Before obtaining an explicit form of P^ 2 \(3;t), we 
prepare the following relations: 

jl 2 \n)\m) 

= ^J(n + l)(n + 2)(m + l)(m + 2)|n + 2)\m + 2), 

(-D - p,A- Ap)\n)\fh) 
= 2nm\n)\m) -2JJn + l)(n + 2)(m + l)(m + 2)|n + 2)|m + 2), 

(A 2 — Op, — \xv) \n) \rh) 
= \J{n + l)(n + 2)(m + l)(m + 2)|n + 2)|m + 2) 

— 2{2nm + n + m + l)|n)|m) + \Jn(n — l)m(m — l)|n — 2)|m — 2), 

(-Z) + z>A + Az>)|n)|m) 
= 2nw|n)|m) — 2\Jn{n — l)m(m — l)|n — 2)|m — 2), 

v 2 \n)\m) 

= \Jn(n — l)m(m — l)\n — 2)|m — 2) for n, m G {0, 1, 2, ...}. (73) 
Using Eq. ( 173"]) . we can write down Pjp(f3;t) as 

1 n/ 2 I ~2 x oq n 

P 2) (M = ^(/3) 2 exp(-^±^)(a|(7|EE 



2 2 ''"^o^oVn! Vm! 

x ([^(n + c) - 2#(n + c + 1) + g(n + c + 2)] 



x^(n + l)(ra + 2)(m + l)(m + 2)|n + 2)|m + 2) 
+ [2nmg(n + c — 1) — 2(2nm + n + m + l)g{n + c) 
+2nmg(n + c + 1)] |n) |m) 

+ c - 2) - 2#(n + c - 1) + #(n + c)] 

x^n(n — l)m(m — l)\n — 2}\m — 2)) 

1 DO OO 

Wexp[-(a 2 + 7 2 )]££ 

Z n=0 m=0 

, a 2n+2 ~2m+2 

x —J—-\g{ Jl + c ) - 2,7(71 + c + 1) + <7(n + c + 2)] 

v nl ml 

a 2n ~2m 

H j ^-2[nmg(n + c + 1) — (2nm + n + m + l)<7(n + c) 

+nmg(n + c — 1)] 

ft 2n-2 r,2m-2 

+ ( w -2)!(m-2)! [g(w + C) " ^ + ^^+9^ + 0- 2)]) 
^(/3) 2 exp(-a 2 )E^r 

Z n=0 n - 

x (Atfa 2 [g(n + c) - 2#(n + c + 1) + g-(n + c + 2)] 
15 
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Figure 2: The graphs of (cr z (t)) given by Eq. (Till and J\{t) given by Eq. ( 1291) for t G [0, 8n] 
and a = 4. (We consider the resonant case, so that we assume Acu = and re — 1.) A thick 
solid curve represents and a thin dashed curve represents (a z (t)). In the numerical 
calculation of (<r z (t)) given by Eq. (JTTJ) , the summation of the index n is carried out up 
to n = 100. In the numerical calculation of J\{t) given by Eq. ( )29|) . the integral / °° dx is 
replaced with J 100 dx. The interval of the numerical integration x G [0, 100] is divided into 
10 6 steps (Ax = 1.0 x 10~ 4 ) and we apply Simpson's rule. To obtain the variation of J\{t) 
against the time t, we divide the interval t G [0, 87r] into 4000 steps (At = 2n x 10~ 3 ) and 
we estimate J\{t) at each time step. Looking at this figure, we notice the following facts. 
In the graph of Ji(t), we can observe only the initial collapse and we cannot observe the 
revival of the Rabi oscillations. [(c z (t)) starts to show the revival of the Rabi oscillations 
around at t = 5tt. By contrast, Ji{t) is nearly equal to zero after t = ir.] 



5 Properties of the integrals that form the atomic 
population inversion 

In this section, numerically evaluating the integrals that form the atomic population 
inversion, we examine their physical meanings. 

First, we consider the case where the electromagnetic field is resonant with the atom at 
zero-temperature, so that we put Au> = and re = 1. We show graphs of J\{t) defined in 
Eq. (T29D and (a z (t)) given by Eq. pi]) in Fig. [2 We show graphs of J 2 (t) defined in Eq. ([29]) 
and (cr z (t)) given by Eq. (flT!) in Fig. [3j In these estimations, we put a = 4. Thus, we can 
neglect the first term of the right-hand side of Eq. ([28]) [—(1/2) exp(— a 2 ) = —5.63 x 10 -8 ]. 

Looking at Fig. El we can conclude that J\ (t) only represents the initial collapse (the 




^(/3) 2 [2(2a 2 7 2 - f - l)P,(t) - 2a\^ + l)Q«(t) 
+4«YQf(t)]. 



(74) 
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Figure 3: The graphs of (cr z (t)) given by Eq. (ITTj) and J 2(f) given by Eq. ( 1291) for t G [0, 87r] 
and a = 4. (We consider the resonant case, so that we assume Atu = and k — 1.) A thick 
solid curve represents and a thin dashed curve represents (a z (t)). In the numerical 
calculation of (a z (t)) given by Eq. (ITTj) . the summation of the index n is carried out up 
to n — 100. In the numerical calculation of J2(t) given by Eq. the integral / °° dy 

is replaced with J 100 dy. The interval of the numerical integration y e [0, 100] is divided 
into 10 6 steps (Ay = 1.0 x 10~ 4 ) and we apply Bode's rule. To obtain the variation of 
J2{t) against the time t, we divide the interval t G [0, 87r] into 4000 steps (At = 2tt x 10~ 3 ) 
and we estimate J2(t) at each time step. Looking at this figure, we notice the following 
facts. In the graph of J2(t), we can observe only the revival and we cannot observe the 
initial collapse of the Rabi oscillations. [(c z {t)) shows the initial collapse from t = until 
around t = n. By contrast, ^(t) is nearly equal to zero from t — to t — 4-7T.] 
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Figure 4: The graphs of (a z (t)) given by Eqs. © and (03JJ) and [1 - 2 exp(-a 2 )/ 1 (0) (t)] 
given by Eq. (1241 for t G [0, 87r], a = 4, Ac; = 4, « = 1 and c = 4. (We consider the 
off-resonant case.) A thick solid curve represents [1 — 2exp(—a 2 )l[°\t)] and a thin dashed 
curve represents (a z (t)). In the numerical calculation of (u z (t)) given by Eqs. (J9]) and (jTOj) . 
the summation of the index n is carried out up to n = 100. In the numerical calculation of 
[1—2 exp(— a 2 )l[°\t)] given by Eq. (I2"4"j) . the integral Jq 00 c&e is replaced with J 100 dx. The 
interval of the numerical integration x G [0, 100] is divided into 10 6 steps (Ax = 1.0 x 10 -4 ) 
and we apply Simpson's rule. To obtain the variation of [1 — 2 exp(— a 2 )l[°\t)] against 
the time t, we divide the interval t G [0, 8tt] into 4000 steps (At = 2tt x 10~ 3 ) and we 
estimate [1 — 2 exp(— a 2 )l[°\t)} at each time step. Looking at this figure, we notice the 
following facts. In the graph of [1 — 2exp(— a 2 )l[ (t)], we can observe only the initial 
collapse and we cannot observe the revival of the Rabi oscillations. [(a z (t)) starts to show 
the revival of the Rabi oscillations around at t — Qir. By contrast, [1 — 2 exp(— a 2 )l[°\t)} 
is nearly equal to Const, given by Eq. (|75|) after t = n.] 

semi-classical limit) and it does not represent the revival (the quantum correction) of the 
Rabi oscillations. By contrast, looking at Fig. [3], we can conclude that only represents 
the revival and it does not represent the initial collapse of the Rabi oscillations. In our 
numerical calculations, we obtain \(a z (t)) — J\(t)\ < 8.0847 x 10 -4 for t G [0,47r] and 
\Mt)) - J 2 (t)\ < 5.6377 x 10- 8 for t G [4^,8^]. 

Here, we pay attention to the representation of Ji(t) given in Eq. ( !29l) . If we let t 
be a large value (t 3> 1), the trigonometric function cos(2y/xt) included in the integrand 
oscillates intensely and rapidly for the small variation of x. Thus, the integral J °° dx of 
J\(t) converges on zero for t ^> 1. Therefore, we can expect that J\(t) never causes the 
revival of the Rabi oscillations. Hence, we can expect that ^(i) lets the revival of the 
Rabi oscillations happen in the range of t 3> 1. 

Second, we consider the case where the electromagnetic field is off-resonant with the 
atom at zero-temperature for estimating the effects of detuning. We show graphs of 
[1 -2exp(-a 2 )/r(t)] defined in Eq. (FJU) and (o z {t)) given by Eqs. © and (JTOJ in Fig.H 
We show graphs of [Const. + 4exp(— a 2 )!^^)] defined in Eq. ( 1241) and (& z (t)) given by 
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Figure 5: The graphs of (cr z (t)) given by Eqs. (jHJ) and ( fTOl) and [Const.+4exp(— a 2 )/^ ^)] 
given by Eqs. §M§ and (J7SD for i G [0, 87r], a = 4, Aw = 4, re = 1 and c = 4. (We consider 
the off-resonant case.) A thick solid curve represents [Const. + 4exp(— a 2 )I^\t)] and 
a thin dashed curve represents (a z (t)). In the numerical calculation of (a z (t)) given by 
Eqs. (Q and ffTUj) . the summation of the index n is carried out up to n = 100. In 
the numerical calculation of [Const. + 4exp( 



-a 



2 ^2°\t)} given by Eqs. (124]) and (1751). 



the integral Jq 00 cfr/ is replaced with / dy. The interval of the numerical integration 
y G [0, 100] is divided into 10 6 steps (Ay = 1.0 x 10~ 4 ) and we apply Bode's rule. To 
obtain the variation of [Const. +4 exp(— at 2 )I^ (t)] against the time t, we divide the interval 
t G [0,8vr] into 4000 steps (At = 2tt x 10~ 3 ) and we estimate [Const. + 4exp(-a 2 )4 0) (t)] 
at each time step. Looking at this figure, we notice the following facts. In the graph of 
[Const. +4exp(-a 2 )4 0) (t)], we can observe only the revival and we cannot observe the 
initial collapse of the Rabi oscillations. [(a z (t)) shows the initial collapse from t — until 
around t = tt. By contrast, [Const. + 4exp(— a 2 
by Eq. (J75D from t = to t = 5vr.] 



j| (t)] is nearly equal to Const, given 
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Figure 6: The graph of P fl M (0; t) given by Eqs. (1711) and (1721) for t G [0, 8n] with putting 
Au = 0, k = 1, c = 0, a = 4, 7 = 1 and 0(/3) = 1/40. 

Eqs. (jUJ) and (fTOl) in Fig. [5], where 



The above Const, is obtained by replacing cos 2 (^x + c\n\t) and sin 2 (\/x + c|/t|t) with 1/2 
in the representation of l{ \t) given by Eq. (j24p . In these estimations, we put a = 4, 
Aw = 4, k = 1 and c = 4. Thus, replacing J °° dx with J 100 (ia;, we obtain Const. = 



Looking at Figs. HJandEl we can conclude that l[ (£) only represents the initial collapse 
and (t) only represents the revival. The reason why (t) given by Eq. ( l24l) only shows 
the initial collapse is as follows. If we let t be a large value (£ 3> 1), we can replace the 
squares of trigonometric functions [cos 2 (\/x + c + l\K\t) and sin 2 (\/a; + c + included 
in the integrand of l[ (t) with their time average 1/2. Thus, we obtain Const, given by 
Eq. ([7S]), and l[°\t) never generates the revival of the Rabi oscillations for t S> 1. 

Third, we consider the resonant case at low temperature. We calculate P^((3;t) 
and Pjj 2 \/3;t) with Aw = 0, so that we estimate the effects of low temperature without 
detuning. Shown in Eqs. fl23|) . ([71]) and (172"]). we can describe the first-order correction of 
the atomic population inversion P^(f3;t) as a sum of l[ (t), li \t), ^(t) and I^it)- 
Similarly, from Eqs. (!23]h ( 172]) and (174]) . we can describe the second-order correction 
Ppifct) as a sum of {l[ l \t) : I E {0, 1,2}} and (4°(t) : ^ <= {0, 1,2}}. From discussions 
given in the previous paragraphs, we can conclude that l[ l \t) generates the initial collapse 
and Iz(t) generates the revival for / = 0,1, 2,.... 

Putting Aco — 0, k = 1, c = 0, a = 4, 7 = 1 and 9{j5) = 1/40, we show graphs of 
Pj} 1 ' (f3;t) and P^ 2 \f3]i) in Figs. |6]and[71 respectively. Looking at Figs. |6]and[71 we find 
that an amplitude of Pjj 2 \[3;t) is comparable to that of P^((3;t). Thus, the parameters, 
7 = 1 and 6((3) = 1/40, are near the boundary of region where the perturbation theory is 




(75) 



-0.2086.... 
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Figure 7: The graph of Pj 2) (/3; t) given by Eqs. and (EI) for t e [0, 8vr] with putting 
Alu = 0, k = 1, c = 0, a = 4, 7 = 1 and 0(/3) = 1/40. 

effective. 

Looking at the graphs of Figs. [6] and [7] and the graphs of Figs. [2] and [3], we understand 
that the effect of P^\(3; t) and Pj 2) (/3; t) on P g (t) is quite small. Hence, we can conclude 
that the collapse and the revival of the Rabi oscillations is robust against the effects of 
low temperature in the region where the second-order perturbation theory is effective. 

6 Discussions 

In this paper, we separate the atomic population inversion of the Jaynes-Cummings model 
into two integrals using the Abel-Plana formula. By numerical calculations, we show that 
the first integral represents the initial collapse (the semi-classical limit) and the second 
integral represents the revival (the quantum correction). Moreover, we examine the time- 
evolution of the JCM with the initial thermal coherent state for the cavity mode at low 
temperature by the second-order perturbation theory. We describe the first and second- 
order corrections as sums of integrals, using the Abel-Plana formula. 

The Abel-Plana formula and its generalized versions are often made use of for the 
calculations of the Casimir energies in different configurations [23j [24j [25] . The author 
thinks that the Abel-Plana formula has a wide application in the field of the quantum 
optics. 
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Figure 8: The closed contour C\ denned on the complex plane. 

A The derivation of the Abel-Plana formula 

In this section, we give details of the derivation of the summation formula described in 
Eq. (|14p . (This formula is introduced in Ref. [IB] without precise derivation. The author 
of Ref. [T5] only mentions a suggestion for proving it.) 

At first, we assume <fi(z) to be analytical and bounded for all complex values of z 
such that n\ < Ke(z) < n 2 , where n\ and n 2 are certain integers. Then, (p(z) / (e~ 2nlz — 
1) is analytical and bounded \/z such that n\ < ~Re(z) < n 2 except for z = n\,n\ + 
l,ni + 2, ...,n 2 — l,n 2 . Hence, if we think about the closed contour C\ shown in Fig. [HJ 
<p(z) / (e~ 2mz — 1) is analytical and bounded inside C\, so that 

-dz = 0. (76) 



Ci e~ 2mz - 1 



(This is the Cauchy integral theorem.) 

Similarly, <p(z)/(e 2mz — 1) is analytical and bounded Wz such that n\ < ~Re(z) < n 2 
except for z = ni, n\ + 1, % + 2, n 2 — 1, n%. Hence, if we think about the closed contour 
C 2 shown in Fig. [91 4>(z)/(e 2mz — 1) is analytical and bounded inside C 2 , so that 

dz = 0. (77) 



fc 2 e 2niz - 1 
Next, we consider the following integrals: 

where C + and C_ are paths on the complex plane shown in Fig. [10J In Fig. [TOl n represents 
an integer such that n\ < n < n 2 and e is a small positive infinitesimal quantity. We can 
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Figure 9: The closed contour C2 defined on the complex plane. 
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Figure 10: The paths C + and C_ defined on the complex plane. Describing C + as an 
opposite path of C+, C + and C_ form a closed contour. 
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Figure 11: The paths Cr + and Cr_ defined on the complex plane, 
rewrite I{n) defined in Eq. (IT8"j) as 



/(B) = _/ 

v ' Jc+ e~ 2mz - 1 Jc 



dz 



c e 2mz - 1 



Jc+ e 2mz - 1 Jc- e 2mz - 1 



* W cfe, (79) 



lc++C- e 2mz — 1 
where C + is an opposite path of C + and $(z) is given by 

_ f e 2 ^0(z) for Im(z) > 
1 > \ cb{z) for Im(z) < ' 1 > 

Inside of the closed contour C + + C_, $(z)/(e 27 ™ 2 — 1) has only one pole, z = n. 
If we write z = n + u, we can expand / (e 27 ™ 2 — 1) around z = n as 

$fn + it) 6(n) 1 „. s . . 

1 ^H, (81) 



e 27ri(n+u) _ J 27TZ W 

where -P(w) is a power series that includes only terms of nonnegative degrees, so that the 
residue is given by 

Res[$(z)/(e 27Tiz - 1), n] = <b(n), (82) 

and we obtain 

J(n) = (b(n). (83) 
Furthermore, we think about the following integrals: 

p-dz+f ^ 



w^L^*! ?£fT*. (84) 



where Cr + and Cr_ are paths on the complex plane shown in Fig. [TTJ In Fig. (TTJ e is a 
small positive infinitesimal quantity. [The index R in Is.(ni), Cr+ and Cr_ implies that 
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Im(z) 

n 2 +ie 




Re(z) 



Figure 12: The paths Cl+ and Cl- defined on the complex plane. 

Cr + and Cr_ are right halves of C + and C_, respectively.] In the limit of e — )■ +0, we 
obtain the following relation: 

Hm Ir(tii) = - lim I{n x ) = -<j){n x ). (85) 
In a similar way, we obtain the following relation: 

, L+ o - - 1 7C L 



IM = L,^i dz+ L ^±i^ (86) 



where Cl+ and Cl- are paths on the complex plane shown in Fig. [12] and 

lim I L (n 2 ) = - lim I(n 2 ) = -</>(n 2 ). (87) 
From Eqs. (175]), (]77]h (]S3]h (^5]) and (|S7]h we obtain 

^ "2— 1 ^ 

-0(ni)+ £ 0(n) + -0(n 2 ) 

Z n=ni+l 

ni+l-e rm+2-e rn 2 -e 



+(/ + / 1 +-+y 1+ ) L-2,L_ 1 + ^fr]^ 

Jni+l+e Jri2 — 1+e C — 1 c — 1 



ni+e 

<-oo ^( n2 _|_ ^) ^ |-n 2 0( x _|_ ^oo) 



A e 2 ^-l y 7ni 



CO 



g— 27ri(a;+joo) ^ 

Looking at Eq. (188]) . we notice the following facts. Because <p(z) is bounded for all z 
such that rti < Re(z) < n 2 , we obtain 

ax = (J, 



"■2 0(a; — ioo) 
m e 2ni ^- ioc ^ - 1 
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cfe = 0. (89) 



The integrand of the fourth term in the left-hand side of Eq. fl88l) can be rewritten as 

<j)(x) 4>{x) _ (e 2nix - l)(j){x) + (e" 2 ™ - l)(j){x) 



Q—2irix ^ ^2-kix ^ ^—2-kix ^^2nix ]\ 

(e 2nix + e~ 2,Tix - 2)(f)(x) 

^ g— 2nix ^Iivix _|_ ^ 

= ~<f>(x), (90) 
so that the limit of the integral, as e approaches zero (e — > +0), is equal to 



rai + 1 — e i>n\+2—e 



n 2- £ <j){x) <f)(x) 



fenW + / + -+/ )[ -11 i + Zix , \ dx 

£^+0 Jni+e Jni+l+e Jn 2 -l+e e zmx — 1 e zmx — 1 

ni 



4>{x)dx. (91) 
Hence, we can rewrite Eq. ( 188|) as 

Z n=m + l z 



4>{x)dx 



hi 

oo 



+f ^ y (n 2 + iy) + + zy) + 0(n 2 - iy) - <ft(m - iy) ^ ^ 



o 



Thus, we obtain the formula of Eq. fTT^j) . If 0(2) — >■ as Re(z) — >■ +00, we obtain 

1 0(0) + £ 0H = r ^(x)^ + i r m ~J { ; iy) dy. 



2 n^l Jo Jo e 2 ^ - 1 

So that, we obtain the formula of Eq. ( TT5l) . This equation is called the Abel-Plana formula. 

B Some remarks about numerical calculations 

In section [5J for the numerical calculations of J\(t), J2{t), Ii(t) and I^it), we use the 
Fortran compiler with quadruple-precision complex (a pair of quadruple-precision real 
numbers). 

We evaluate the Gamma function included in Eqs. ff24|) and (1291) numerically by the 
Lanczos approximation 



\nT{z) = (z + ^)\n(z + g + ^)-(z + g + ^) 



+ ln[^(c + V^^ + e)] forRe(z)>0, (94) 



where g = 5, 

c = 1.000 000 000 190 015, 
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ci = 76.180 091 729 471 46, 

c 2 = -86.505 320 329 416 77, 

c 3 = 24.014 098 240 830 91, 

c 4 = -1.231 739 572 450 155, 

c 5 = 0.120 865 097 386 617 9 x 10~ 2 

c 6 = -0.539 523 938 495 3 x 10~ 5 . (95) 

The approximation with Eq. (I94"|) gives the error upper bound |e| < 2 x 10 -10 . 

For carrying out the numerical integration of J\{t) defined in Eq. ( 12 9 p and I\\t) 
defined in Eq. ([241 . we use Simpson's rule [26]. For carrying out the numerical integration 
of J 2 (t) defined in Eq. ([29]) and 4°(t) defined in Eq. we use Bode's rule [26]. 

In the numerical integration of J2{t), we pay attention to the following fact. In the 
limit as y — > +0, the integrand of converges on a finite value as 

1 \a\ 2iy 1 — 1 

]%> ^^l lm{ rTTT^) C °< 2 ^ = 2^(21b H + 7 - 2t 2 ), (96) 

where 7 = 0.577215... is the Euler-Mascheroni constant. Similarly, in the limit as y — > +0, 
the integrand of (t) converges on a finite value as 

2iy 



1- — -lm{— ——[cos 2 (Jc + iy\K\t) + — — sin 2 ( Jc + iy\K\t)}} 
y->+o e Z7ry — 1 l{l + ty) v c + iy v 

= - — [-l + 2c7 + cos(2Vck|t) +4cln|a|]. (97) 

4C7T 

As shown above, when we calculate I 2 (t) numerically, we have to be careful about taking 
the limit of the integrand of I-Pfy) as y — > +0. 

The numerical integration of ^{t) is more difficult than that of J\{t). The numerical 
evaluation of J2{t) for t > 8rr does not converge on a reasonable value even if we use 
Romberg's method [26]. The exactly same things happen when we calculate l[°\t) and 
l2°\t) numerically. 
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